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ABSTRACT 

Observations are presented of emission lines from organic molecules at frequencies 32 - 50 GHz in the 
vicinity of Chamaeleon MMS1. This chemically -rich dense cloud core habours an extremely young, very 
low-luminosity protostellar object and is a candidate first hydrostatic core. Column densities are derived and 
emission maps are presented for species including polyynes, cyanopolyynes, sulphuretted carbon-chains and 
methanol. The polyyne emission peak lies about 5000 AU from the protostar, whereas methanol peaks about 
15,000 AU away. Averaged over the telescope beam, the molecular hydrogen number density is calculated to 
be 10 6 cirT 3 and the gas kinetic temperature is in the range 5 - 7 K. The abundances of long carbon chains 
are very large, and are indicative of a non-equilibrium carbon chemistry; and HC7N column densities 
are 5.9^"| x 10 n cirT 2 and 3.3^j ( j x 10 12 cirT 2 , respectively, which are similar to the values found in the 
most carbon-chain-rich protostars and prestellar cores known, and are unusually large for star-forming gas. 
Column density upper limits were obtained for the carbon-chain anions C4H and C(,H , with anion-to-neutral 
ratios [C4H~]/[C4H] < 0.02% and [C6H~]/[CeH] < 10%, consistent with previous observations in interstellar 
clouds and low-mass protostars. Deuterated HC3N and 0C3H2 were detected. The [DC3N]/[HC3N] ratio of 
approximately 4% is consistent with the value typically found in cold interstellar gas. 

Subject headings: Astrochemistry - Stars: formation - ISM: abundances - Radio lines: ISM, stars 



1. INTRODUCTION 

During the earliest stages of star formation, quiescent 
molecular gas collapses to form a dense, prestellar core and 
subsequently a protostar with accretion disk, from which 
planetary systems may eventually arise. Understanding the 
evolution of molecular complexity as matter passes through 
these phases is crucial to the development of theories concern- 
ing organic chemistry during the formation of planeta ry sys- 
tems, with relevance to the origin of life in the Galaxy (Ehren- 
freund & Charnley 2000). Interstellar clouds and prestellar 
cores have been shown to host a rich organic chemistry that 
results in the synthesis of a ran ge of large organic molecules 
( |Herbst &~v an Dishoeck 2009). Unsaturated carbon-chain- 
bearing species such as the polyynes and cyanopolyynes are 
high ly abundant in dark inte rstellar clouds such as TMC-1 
{e.g. Cernicharo et al. 1984) ), and have recently been found 
to be present in large concentrat ions around some prestel- 
lar cores and low-mass protostars ( Sakai et al.|200 8b 2009a; 



Gupta et al.|2009| |Cordiner et al.|201 \) . It has been hypoth- 
esised that such complex molecules may survive the passage 
from protostellar envelope to accretion disk, to later become 
incorporated, largely un modified, into solar system bodies 
such as comets (e.g. [Meier et al.|199 8 ; Willa cyl2007| |Charn-| 
ley & Rodgers|2008||20097 lMumma & Char nley|201 A) . Fur- 



ther observational and theoretical studies of the chemical evo- 
lution of gas and dust in the vicinity of protostars as they 
form and heat their surroundings will be crucial in order to 
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understand the chemical reagents present at the start of planet 
formation, and those delivered to planetary surfaces later by 
cometary impacts. 

Chamaeleon MMS1 (abbreviated Cha-MMSl) is an ideal 
source in which to study the chemistry of a protostar in its 
natal environment. Thus far, Cha-MMS 1 has been relatively 
little-studied in terms of its molecular inventory. It was first 
detected in 1.3 mm emission as a ve ry compact object (=s 30" 
in diameter) by Reipurth et al. ( 1996 1. Infrared and sub-mm 
photometry are consistent with the presence of a young Class 
protostar embedded in a dusty envelope with a core den- 
sity of 3.3xl0 6 c m" 3 and diameter of about 3600 AU ( Bel- 
lloche et al.||20061). The central object has been classified by 



1120061 

|Belloche et al.| ( |2006 |201 \\ as either an extremely young, 



very low-luminosity Class protostar (VeLLO) or an exam- 
ple of a first hydrostatic core (FHSC). VeLLOs are likely to be 
young protostars undergoing episodic accretion, and a number 
of these were identified by the recent Spitzer c2d survey ( |Dun-| 
|ham et al.|20 08j). The FHSC is a transi ent, quasi-hyd rostatic 
object predicted by the simulations of Larson ( 1969), but so 
far eluding observational confirmation. It is theorised to be 
formed in the early stages of protostellar collapse at the point 
when the density becomes sufficiently high that the core is 
opaque to thermal radiation. Virial heating then provides the 
energy to pressurise the (predominantly molecular hydrogen) 
gas against further collapse. Representative of a phase of pro- 
tostar evolution preceding Class , the FHSC is also known 
as Class -I (Boss & Yorke 1995 1. The physics and chem- 
istry of FHSCs and VeLLOs have recently been discussed by 



Omukai ( 2007 1 and Lee ( 2007 ), respectively. The predicted 
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Fig. 1 . — Example 7 mm Mopra spectra of Cha-MMS 1 . For CgH , the 
spectrum shown is the average of the three observed lines. 

spectral energy distribution of an FHSC is similar to that of 
the coolest, lowest-luminosity VeLLOs, and at present there 
are only four likely candidate FHS Cs reported in the liter- 
ature: C ha-MMS 1, L1448 I RS2E ( |Chen et al.||20l0), Per- 
Bolo 58 ( |Enoch et al71[20TOl > and LI 451 -mm ( [Pineda et aL] 
|20TT) . 

Observations of organic molecules in Cha-MMS 1 were 
made previously b y |Kontinen et al.| ( |2000| ) using SEST; Bel- 
|loche et aL] ( |2006| l observed a selection of deuterated and un- 
deuterated molecular ions, and |Tennekes et alT| ( |2006[ ) mapped 
the HCN and HNC emission around the protostar. However, 
the chemical and physical characterisation of Cha-MMS 1 is 
still far from complete. As part of a program to understand 
organic chemistry during the earliest phases of low-mass star 
formation, we have initiated molecular line surveys in a sam- 
ple of around 70 young protostars, focussing on individual 
sources for more detailed followup. In this article, we present 
7 mm spectral line data and maps of cold carbon-chain- 
bearing species and other organic molecules in the vicinity 
of Cha-MMS 1. These observations are analysed in the con- 
text of the physical and chemical properties of the protostellar 
envelope which, compared with other young carbon-chain- 
rich protostars (such as L1527 and L1521F, whose o rganic 
chemistry were recen tly examined by Sakai et al. 2008b and 
Takak uwa et al.|201 1[ ), help to constrain the current state and 
the history of the core. 

2. OBSERVATIONS 

Observations of Cha-MMS 1 were carried out in March 
2010 using the ATNF Mopra 22 m radio telescope. The Mo- 
pra Spectrometer (MOPS) was configured in zoom mode to 
simultaneously observe 16 spectral windows distributed over 
an 8 GHz range, each with a bandwidth of 137.5 MHz and a 
channel spacing of 33.6 kHz. During observations, the total 
system temperatures were in the range 80 - 200 K. The beam 
efficiency factor at 7 mm (42 GHz) was 0.6. Pointing was 
checked about every three hours and was typically accurate to 
within 5 - 10". 

Position-swit ched integrations of Cha-MMS 1 were per- 
formed at the ( jKontinen et al.|2000) > peak HC3N position, at 
(J2000) coordinates R.A. = 1 1:06:33, deck = -77:23:46. The 
off position used was R.A. = 1 1:17:38, decl. = -77:23:25. To 
help mitigate the known Mopra 30 MHz spectral baseline rip- 



ple, the telescope subreflector was continuously driven up and 
down during observations, with a displacement of 3.5 mm. 
Atmospheric opacity corrections e T " are required at 7 mm; 
zenith optical depths as a function of frequency (tq(v)) were 
calculated using the GILDAS ATM routina^lbased on the pre- 
cipitable water vapour column and atmospheric temperature 
and pressure at the time of observations. 

Observations were performed using two different MOPS 
configurations centered at frequencies of 36 GHz and 45 GHz, 
respectively. Approximately 22 hr of on-source integration 
was obtained for the 36 GHz setting, yielding RMS main- 
beam brightness temperature sensitivities 7' MB (RMS) ~ 6 mK 
per channel. Two hours integration in the 45 GHz setting 
yielded Tmb(RMS) ~ 30 mK per channel. The observed 
molecular transition lines are given in Table [T] The telescope 
beam FWHM was 96" at 36 GHz and 77" at"45 GHz. 

On-the-fly mapping of the 8' x 8' region surrounding Cha- 
MMS 1 was performed using the 36 GHz setting, employing 
a Nyquist sampling rate in the R.A. (scan) direction and 33" 
spacing between map rows. The resulting map sensitivity was 
I'mb(RMS) ~ 70 mK per channel per pixel. 

Spectral line data were reduced and combined using stan- 
dard routines in the ATNF Spectral Analysis Package (ASAP) 
version 3.0. Baseline subtraction and spectral line analysis 
were performed using NRAO IRAFF] Examples of a selection 
of the observed spectral lines are shown in Figure [T] Map 
data were reduced usi ng the AIPS++ p rograms LIVEDATA 
and GRIDZILLA (see Barnes et al. 2001), and smoothed with 
a Gaussian FWHM of 96" applicable to the beam size at 
36 GHz. 

3. RESULTS 

3.1. Molecular emission maps 

Maps of the integrated line intensity J dv are shown for 
HC 3 N, HC 5 N, C 4 H, CH3OH, CCS and C 3 S in Figures [2] to 
[4j All mapped species show resolved structure with strong 
emission peaks within 2' of the protostar (which is plotted as 
an asterisk). The polyyne-based species HC3N, HC5N and 
C4H all show peaks that are approximately spatially coinci- 
dent with each other and are elongated in the east-west direc- 
tion. The sulphuretted carbon chains show a similar morphol- 
ogy, but are offset by about 35" to the west of the polyynes. 
Elliptical Gaussian fits to these carbon-chain-bearing species 
have FWHM ~ 200" x 400", with the long axis orientated at 
a mean position angle -5° from the horizontal (E-W) axis. 
The methanol peak has similar dimensions, but is located 
70" south-west of the polyyne peak. Extended structures are 
present in all the carbon chain maps, the principal common 
feature being a ridge in the east-west direction, the full extent 
of which cannot be determined from our maps. The methanol 
distribution is markedly different however, with relatively lit- 
tle extension to the east of the protostar, and a secondary 
peak to the north-east. The polyyne-based species peak clos- 
est to the protostar (about 30 - 40" away), followed by the 
sulphuretted carbon chains (60" away), then methanol (100" 
away). Assu ming a distance to Cha-MMS 1 of 150 pc (Knude 
& Hog|[T998| l, these correspond to respective distances from 
the protostar of about 5,000 AU, 9,000 AU and 15,000 AU 

2 http://www.iram.fr/IRAMFR/GILDAS 

3 IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 
omy (AURA) under cooperative agreement with the National Science Foun- 
dation. 



7 mm spectroscopy of the young protostar Chamaeleon MMS 1 
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TABLE 1 
Observed line parameters 



Species 


Transition 


Freq. 




VLSR 


Av 


/ 7/ M b dv 






(MHz) 


(mK) 


(kms~') 


(kms^ 1 ) 


(mKkm s _1 ) 


C 6 H 


2 n 3/2 7 = 23/2 - 21/2 / 


31881.8596 


17(7) 


4.71 (15) 


0.93 (34) 


17(6) 


C 6 H 


2 n 3/2 7 = 23/2- 21/2 e 


31885.5414 


21 (6) 


4.58 (13) 


1.08 (31) 


25 (7) 


C 3 H 


7 = 3/2- 1/2 eF = 2- 1 


32660.6550 


86 (10) 


4.75 (2) 


0.78 (6) 


72 (8) 


C 3 H 


7 = 3/2- 1/2 eF = 1 -0 


32663.3750 


39(7) 


4.86 (5) 


0.79 (14) 


32(6) 


C 3 H 


7 = 3/2- 1/2 eF = 1 - 1 


32667.6366 


22(7) 


4.37 (12) 


0.86 (27) 


20 (6) 


HC 7 N 


7 = 29 - 28 


32711.6820 


19(11) 


4.70 (9) 


0.47 (24) 


14(6) 


C 6 FT 


7 = 12- 11 


33044.4879 


<16 






<13 


CC 34 S 


Nj = 2 3 - 1 2 


33111.8390 


34(7) 


4.76 (7) 


0.94 (16) 


37 (7) 


CCS 


JV> = 2 3 - 1 2 


33751.3737 


792 (45) 


4.81 (1) 


0.76(1) 


650 (38) 


DC 3 N 


7 = 4-3 


33772.5312 


58(7) 


4.71 (5) 


1.3(10) 


80 (9) 


HC 7 N 


7 = 30 - 29 


33839.6250 


<19 






<16 


HC 5 N 


7=13-12 


34614.3870 


368 (22) 


4.76(1) 


0.77 (1) 


301 (19) 


C„H 


2 n 3/2 7 = 25/2 - 23/2 / 


34654.0373 


13(6) 


4.68 (22) 


1.28(48) 


22(6) 


C 6 H 


2 n 3/2 7 = 25/2- 23/2 e 


34658.3831 


14(5) 


4.71 (20) 


1.39 (49) 


18(5) 


C 3 S 


J = 6-5 


34684.3675 


247 (17) 


4.70(1) 


0.76 (2) 


200 (15) 


C 6 H 


2 n l/2 J = 27/2 -25/2 / 


34887.1145 


<18 






<15 


C 6 H 


2 ni /2 / = 27/2- 25/2 e 


34917.8847 


<18 






<15 


HC 7 N 


J = 31 - 30 


34967.5880 


11 (6) 


4.75 (27) 


1.52 (65) 


17 (5) 


C„H- 


7 = 13 - 12 


35798.1532 


<17 ■■■ 






<14 ■■■ 


C 4 H 2 


Jk„k c = 4i3 - 3 12 (ortho) 


35875.7746 


33 (6) 


4.65 (7) 


1.14(18) 


41 (7) 


HC 7 N 


7 = 32-31 


36095.5340 


24 (12) 


4.67 (7) 


0.39 (17) 


10(4) 


CH 3 OH 


7^ = 4_] - 3 E 


36169.2610 


94 (10) 


4.71 (2) 


0.87 (5) 


86 (9) 


HC 3 N 


7=4-3F=4-4 


36390.8861 


76 (9) 


4.77 (3) 


0.82 (6) 


83 (9) 


HC 3 N a 


7 = 4-3AF = -l 


36392.3240 


1673 (90) 






2158 (114) 


HC 3 N 


7 = 4- 3F = 3-43 


36394.1765 


83 (9) 


4.74 (2) 


0.74 (5) 


65 (7) 


HC 7 N 


7 = 33-32 


37223.4920 


24(7) 


4.82 (9) 


0.67 (20) 


19 (6) 


C 4 FT 


7 = 4-3 


37239.4020 


<17 






<14 


HC 5 N 


7= 14-13 


37276.9940 


264(18) 


4.77(1) 


0.79 (2) 


223 (16) 


C6H 


2 Tl-xn J = 2112 - 25/2 f 


37426.1916 


31 (8) 


4.87 (6) 


0.69 (15) 


24 C51 




2tTo ^ f — 77 /9 — 7S /7 ^ 


1,7 A Q1 7SS0 
J /*+3 1 .ZJJU 


1 Q (^\ 

iy (3) 


A f\ 1 ( 1 Ai 
4.01 yi'r) 


1 S (11} 
1 .j \ J3 } 


39 (1\ 


C4H 


N — 4—3 / — 9/9 — 7/7 


38049 6583 


398 (24) 


A 8Q (1 ) 


1 1 5 n 1 

L.L3 1,1,1 


494 (32) 


p.u 
^4^1 


N — A 3 /— 7/7 — ^/7 


38088 4S87 


38S (931 


/ \ l } 


n cn n ~i 

V.y3 \1 } 


38/1 (1A\ 
J04 iz-^t ) 


ru n run 


J K a K c — 2(12 — 2oi E 


38SOS QQQft 


AA ifi\ 


A SO 1S1 


n G7 fi n 
u.o / 1117 


Aft (f^t 


fu. run 


^KaKc — 2o2 ~ 2oi A 
/ — 1 4 11 


38 VI 7 1 130 
JOJ li.l 1JU 


SQ 181 


A Q7 /31 
4.y/ (3) 




33 (1) 




38ssi 8083 


^17 






<14 


nr.w 

U\^ 3 1N 

HCO + 


7 — S — A 


4971 S S877 


1 08 (\ 71 




v.yj 11Z 7 


1 1 n n si 

11U 1 1 J 


JK a K c — 2q2 — lot 


A11{*(\ 1 Q7S 


^Sl 






<42 


* / ■ 1 1 ('11 ( \ 1 1 


7 — 1 ft 

-/s: a: c - 'oi - Ooo 


/i ift7£ /^ftnft 


<J 1 






7 

<4Z 


DC 5 N 


7= 17- 16 


43215.4382 


<51 






<42 


HNCO 


JlC a K c — 2q2 — lot 


43963.0000 


100 (24) 


4.45 (7) 


0.72 (17) 


81 (24) 


CCS 


N, = 4 3 - 3 2 


43981.0270 


127 (30) 


4.70 (5) 


0.54(11) 


86 (20) 


c-C 3 H 2 


Jk q k c = 321 - 3i2 (ortho) 
N = 5 -4 7= 11/2-9/2 


A A 1 f\A "IOl~\£. 

44104. /806 


189 (30) 


4.48 (3) 


0.63 (8) 


128 (21) 


C 4 D 


44146.0772 


<63 






<52 


C 4 D 


/V = 5-47 = 9/2-7/2 


44182.0627 


<63 






<52 


HC 5 N 


7= 17- 16 


45264.7199 


278 (34) 


4.53 (2) 


0.62 (5) 


193 (25) 


HC 3 N a 


7 = 5-4AF = -1 


45490.3138 


2198(132) 






1965 (115) 


c-C 3 H 2 


/x„x c = 2n - 2 2 (para) 


46755.6136 


840 (63) 


4.56(1) 


0.64 (2) 


618 (54) 


r-CH 3 CH 2 OH 


■/iCA = 4o4 - 3l 3 


46832.8000 


<81 






<67 


C 4 H 


/V = 5-47 = 11/2-9/2 


47566.8090 


635 (54) 


4.68(1) 


0.64 (3) 


443 (38) 


C 4 H 


N = 5-47 = 9/2-7/2 


47605.4850 


515 (46) 


4.55 (1) 


0.65 (3) 


380 (39) 


HC 5 N 


7= 18- 17 


47927.2746 


226 (43) 


4.52 (3) 


0.45 (7) 


106 (21) 


CH 3 OH 


7* = In - A + 


48372.4600 


657 (68) 


4.51(1) 


0.52 (3) 


365 (37) 


cs 


7=1-0 


48990.9549 


2465 (145) 


4.44(1) 


0.91 (1) 


2579 (160) 


C 3 N a 


iV = 5-47 = 11/2-9/2 


49466.4200 


117 (42) 


4.72 (9) 


0.68 (24) 


77 (21) 


C 3 N 


N = 5-47 = 9/2-7/2 


49485.2480 


<126 






<104 


c-C 3 HD 


^/v,Ai = hi -Oqo 


49615.8517 


366 (54) 


4.48 (3) 


0.56 (6) 


222 (34) 


Note. — Errors on the last quoted digit(s) are given 


in parentheses. 


Upper limits are 


3Ai> X RMS noise, with Av 


= 0.78. 



d HC3N and C3N measurements are integrated over the three main hyperfine components. 
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in the plane of the sky, and indicate a gradient in the chemi- 
cal/physical properties of the gas. 

3.2. Molecular column densities 



TABLE 2 

Molecular column densities and abundances 



All spectral lines - observed at the Kontinen et al.| 
HC3N peak - were fitted using single Gaussians. Integrated 
line intensities and fit parameters are given in Table [TJ One- 
sigma errors were calculated using Monte Carlo noise re- 
sampling with 1000 replications. The errors in Tmb and 
J Tmb dv given in the table also include uncertainties of +5% 
to account for errors in T MS introduced by the atmospheric 
opacity correction as a result of uncertainties in the atmo- 
spheric model parameters. Measured line FWHM are in the 
range 0.4 - 1.5 kms -1 . Unresolved hyperfine structure is 
present in many of the lines, and causes significant addi- 
tional broadening of the observed line widths for HC3N, C3N, 
DC3N, C4H and CeH. Ignoring outliers more than one stan- 
dard deviation from the mean, the average line FWHM is 0.78 
kms -1 . 

Molecular column densities are presented in Table|2] These 
were calculated using a variety of methods as indicated in the 
table notes. Spectral line data and partition functions have 
been taken from The Colo gne Database for Molecular Spec- 
troscopy (Miiller et al.pOOll , with supplementary data from 
the JPL databasfl Partition functions for the different nuclear 
spin configurations of 0C3H2 and CH3OH were calculated by 
direct summation over the relevant energy levels. The beam- 
filling factor was assumed to be unity for all emission, which 
is justified based on the spatial extent of the emission maps. 
However, structure on size scales less than the size of the Mo- 
pra beam is likely given the presence of the compact proto- 
stellar core. This may introduce additional uncertainties into 
the derived column densities and temperatures. 

3.2.1. Rotational diagrams 

Rotational diagrams (see, e.g. Cummi ns et al.| [l986) were 
used for those species for which at least two transitions from 
significantly different upper-state energies were observed. 
The rotational diagrams containing at least three transitions 
each are shown in Figure [5] Where possible, in order to bet- 
ter constrain the temperatures and column densities, our data 
were supplemented w ith data from the SEST observations of 
Kontinen et al. (|2000[). Derived rotational excitation temper- 



atures T rot are relatively low (< 10 K), so it was necessary to 
account for the cosmic microwave background (CMB) radia- 
tion by including a factor [1 - (J v (Ti, s )/ J v (T rot )]~ l in the col- 
umn density calculations, where J Y (T) is the Planck radiation 
law and T bg = 2.73 K is the CMB temperature. For HC 3 N, 
C4H, ortho-c-C3H2 and CCS, the peak opacities (r v ) are sig- 
nificant for one or more lines, so for these species the points 
on the rotational diagram were corrected using the factor 

Cv = T v /(l-e- T '), (1) 
where r v was derived from the opacity equation 

^MB = UviTrot) ~ Jy(T bg W ~ e~ T ") (2) 

The correction factors were determined iteratively, with ini- 
tial values (c°) calculated at temperatures of T® ol (HCi,N) = 
7.6 K, r r o( (C 4 H) = 4.3 K, r° ,(ortho-c-C 3 H2) = 7.7 K and 
r ; o ,(CCS) = 5.8 K, which were derived from the respective 

4 http://spec.jpl.nasa.gov/ 



Species 


N (cm -2 ) 


T IK) 


Abundance 


Method 11 


C 3 H 


2.9+0;' (12) 


6.1 ± 1.4 


9.5+i° 6 (-ll) 


(1) 


C 4 H 


1.4+gJ(14) 


4.9 ± 1.0 


4.6+|« (-9) 


(2) 


C4H" 


< 4.5(10) 


6.1 ± 1.4 


< 2.4 (-12) 


(1) 


C 6 H 


5.9+ff (11) 


6.1 ± 1.4 


2.0+ 2 J (-11) 


(1) 


C 6 H~ 


< 4.8 (10) 


6.1 ± 1.4 


< 2.5 (-12) 


(1) 


C3N 


2.5+^(12) 


6.1 ± 1.4 


8-2+' 3 3 (-11) 


(1) 




z.U {Lj) 




O. / \— 1U) 


(i) 




1 a+0.2 /1 t\ 
1.0 _ , (13) 


/ .z ± u.z 


^ -3.8 (-1U ' 


I A \ 




A C+1.6 / ] \ 

4 -1.2 


"7 1 j. c\ s 
/ . 1 ± U.o 


1 c+ 1.7 / 1 n\ 


(z) 




a -2+8.O /1 n\ 


0. 1 ± 1 .4 


1 1 +4.9 / 1 a\ 
lA -QJ 


/ 1 

( i) 


DC3N 


"7 "3+2.3 /I 1 \ 

'■ 3 _l.l C 11 ) 


6.1 ± 1.4 


a +2.6 / 1 1 \ 

2-4_ 9 (-11) 


(2) 




< 4.V (,12) 


D.l ± 1.4 


< Z.O \— L\J) 


(1) 


O-C-L-3 rn 


1 q+5.9 1 1 t\ 
1,J -0.4 


A 1 -1- 1 A 
0. 1 ± 1 .4 


a o+31 t \r\\ 

4 - J -2.o y~ l[J ) 


iT\ 
\i) 


u-c -L.3112 


7+O.9 /1 9-1 


f. f. _i_ A O 

0.0 ± U.J 


9 Q+2.1 / 1 n\ 


(A\ 
W 


P-C-C3H2 


1 A + l.O / 1 1\ 

l-4_ 03 (13) 


6.1 ± 1.4 


A T + 7.3 / 1 f\\ 

4.7 2 (-10) 


(3) 


P-0C3H2 


1.5+1;' (13) 


6.1 ± 1.4 


5.0+5-5 ( _ 10) 


(1) 


0C3HD 


4.5+i;3 (12) 


6.1 ± 1.4 


1.5^5 (-10) 


(1) 


C 4 H 2 b 


2.1+3j(l2) 


6.1 ± 1.4 


6.8! 23 „(-ll) 


(1) 


HNCO 


7.6+IJ (12) 


6.1 ± 1.4 


2.5+ 3 ; 3 (-10) 


(1) 




< 9.0 (1 1) 


6.1 ± 1.4 


< 4.7 (—1 1) 


(1) 


CH 3 OH-A 


6.0+19(12) 


7.1 ±0.9 


2.0i 2 i(-10) 


(4) 


CH30H-E 


1.2+g | (13) 


4.7 ±0.6 


4.0+|3 (-10) 


(4) 


CH3CHO 


3.6+10(12) 


6.1 ± 1.4 


1.2+1; 2 (-10) 


(1) 


(-CH 3 CH 2 OH 


< 1.2(13) 


6.1 ± 1.4 


< 5.0(-10) 


(1) 


CS 


> 3.8(13) 


6.1 ± 1.4 


> 9.5 (-11) 


(1) 


CCS 


1.3+° : 2 (13) 


5.6 ±0.2 


4.2+j 3 (-10) 


(4) 


cc 34 s 


7.5+^(11) 


6.1 ± 1.4 


2.5!^ (-11) 


(1) 


C 3 S 


2.3+0 9(12) 


6.1 ± 1.4 


7.6^ (-11) 


(1) 



Note. — Base-ten exponents are given in parentheses. Col- 
umn density upper limits are 3<x. Abundances are given with 
respect to H 2 column density (3 ± 1) x 10 22 cm 4 . 

a Method used for column density calculation: (1) LTE with 
T rot = 6. 1 ± 1 .4; (2) Rotational diagram; (3) Best-fit RADEX ra- 
diative transfer model (with density set to nn 2 = 10 s cm -3 ); (4) 
Rotational diagram including supplementary data from ( Konti- 
|nen et al.|2000| . 

Total C4H2 column density calculated from ortho species as- 
suming an equilibrium ortho-to-para ratio (3:1). 

uncorrected rotational diagrams for these species. For each 
successive iteration, new correction factors (c") were calcu- 
lated based on the temperatures derived from the previous it- 
eration. Converged rotational temperatures and correspond- 
ing column densities are given in Table [2] 

3.2.2. LTE calculations 

For molecules for which rotational diagrams could not be 
plotted, column densities were calculated under the assump- 
tion of local thermodynamic equilibrium (LTE), using a gas 
temperature of 6.1 + 1.4 K, which is the mean average (±lcr), 
of the excitation temperatures derived from the rotational dia- 
grams. 

For molecules in this category with multiple line detections, 
the weighted-average column densities were calculated based 
on the individual line signal-to-noise ratios. Where molec- 
ular lines were searched for but not detected, column density 
upper limits were calculated using upper integrated line inten- 
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Fig. 4. — Integrated emission intensity maps of CCS Nj = 4^ - 32 (left) and C3S 7 = 6-5 (right). The protostar Cha-MMSl is marked with an asterisk, which 
denotes the center of the millimeter source. Telescope beam shown in bottom left. 



6 



Cordiner et al. 



sity limits of 3Av x (RMS noise), for which Av = 0.78 km s" 1 
was used. For C(,H , the upper limit of 4.79 x 10 10 crrT 2 was 
derived from the average of the three observed C6H spectra. 

3.2.3. Radiative transfer modelling 

For HC3N and 0C3H2, the observed emission lines were 
subject to radiative transfer modelling using the RADEX code 



ing to greater line opacities a nd column densities. For ex- 
ample, |Kontinen efaT] p000~| derived T ex = 4.7 ± 0.2 K, 



developed by van der Tak et al. (2007). This routine employs 



statistical equilibrium calculation for molecular excitation in- 
volving collisional and radiative processes, and accounts for 
optical depth effects using an escape probability method. Col- 
lisional and radiative (de-)excitation rates were taken from the 
Leiden Ato mic and Molecular Database (LAMDA) ( |Schoier| 
et al.| |2005 P] which tabulates scaled versions of the original 



data pu blished by |Green & Chapman (1978 ) for HC3N and 
|Chandra & Kegel] ( |2000| l for c-C 3 H 2 . The RADEX model 
free parameters (number density of primary collision partner 
«h 2 1 g as kinetic temperature T and molecular column density 
N) were optimised using a least-squares algorithm to produce 
the best fit to the observed integrated line intensitie s (includ- 
ing the add itional higher- frequency line data from |Kontinen| 
let al.|2000] >. 

The two nuclear-spin species of 0C3H2 (ortho and para) 
were considered separately. For ortho-c-C3H2, the following 
transitions were included in the fitting: Jk„k c = 2i2~ loi> 3i2- 
221 and 32i - 3i2- The best-fitting density was 1.0 x 10 6 cm -3 , 
but T and N were not well constrained by the data. Therefore, 
in addition to the LTE values, RADEX column densities for 
ortho- and para-c-C3H2 are given in Table [2] calculated us- 
ing a fixed density of n Hl = 1.0 x 10 6 cirT 3 and temperature 
r = 6.1±1.4K. Even at densities ~ 1 6 , the complexity of the 
C-C3H2 spectrum results in significant departures from LTE 
excitation, so the results from radiative transfer modelling 
are preferred over the LTE /rotational diagram results. How- 
ever, the calculated ortho-c-C3H2 energy level populations are 
highly dependant on the assumed temperature, which leads to 
substantial uncertainty in the column density. On the other 
hand, the HC3N rotational excitation is close to LTE, as shown 
by the near-linearity of the corrected rotational diagram in 
Figure [5] As a result, the gas density is only constrained by a 
lower limit («# 2 > 5.5 x 10 5 cm -3 ), derived using a RADEX 
fit to the J = 4-3, 5-4, 10-9 and 11 - 10 transitions. Above 
this value, the best-fitting HC3N column density is relatively 
density-insensitive; the value for A^HQN) shown in Table [5] 
was calculated for n# 2 = 1.0 x 10 6 crrT 3 . The best-fitting 
HC3N excitation temperature was 7.1 K, which matches well 
the value of 7.2 + 0.2 K derived from the corrected HC3N 
rotational diagram. 

4. DISCUSSION 
4.1. Molecular column densities and abundances 



Based on the dust emission, Tenneke s et al.| (|2006) calcu- 
lated an average H2 column density of N{H2) = (3 ± 1) X 
10 22 crrT 2 in a 50" beam surrounding Cha-MMSl. Using this 
range, the observed molecular abundances relative to H2 are 
given in Table [2] 

Our column densities for CC S c-C 3 H 2 , HC 3 N and C H 3 OH 
differ from those published by Kontinen et al. ( |2000| l, who 
used a line-ratio method for the derivation or LTE - excita- 
tion temperatures. Their derived temperatures are signifi- 
cantly less than ours (for all species except A-methanol), lead- 

5 http://www.strw.leidenuniv.nl/~moldata 



N = 4 7 X 10 13 cm" 2 forT^C?^, and T ex = 4.1 ± 0.1 K, 
N - 4.5 x 10 14 cnr 2 for HC3N. For their strongest observed 
lines of these species, core opacities of r v w 3 were derived, 
whereas we derive opacities less than 1 based on their data 
for the same transitions. Their HC3N column density, in par- 
ticular, seems unrealistically large compared to other carbon- 
chain rich interstellar clouds (even the TMC-1 cyanpolyyne 
peak); we derive values of 2.0 x 10 13 cm -2 using RADEX 
modelling and 1.6 x 10 13 cm" 2 using a (corrected) rotational 
diagram. Their use of approximate partition functions for or- 
tho and para 0C3H2 is another source of discrepancy. 

The mean gas number density within a FWHM = 50" 
beam surrounding Cha-MMSl was found by Belloche et al. 
( 201 1[ ) to be 9.8 x 10 5 crrT 3 (assuming a dust temperature of 
2 K), which closely matches our derived H2 number density 
of riH 2 = 1.0 x 10 6 cm" 3 . If the dust temperature is closer 
to our mean derived value of 6.1 K, the density calculated 



by Belloche et al.| ( 201 1 1 could be up to about a factor of 3 
greater. This would still be consistent with our results how- 
ever, given their smaller beam size, which probed more of the 
higher-density material nearer to the protostar. 



4.1.1. Polyynes and cyanopolyynes 

The abundances of carbon-chain-bearing species (including 
the polyynes and cyanopolyynes) in Cha-MMSl are large, 
and similar to those found in other carbon-chain-rich inter- 
stellar clou ds, prestellar cores and protostars suc h as TMC - 
1, L1527 dSakai et al.||2008b"ll, L1521E (|Hirota et al.||2004l ), 
L1512 and L1251A ( [Cordiner et al.|201 1) . 

Previous laboratory and theoretical studies have shown that 
polyynes and cyanopolyynes may be rapidly synthesised in 
cold, dark interstellar clouds throug h gas-phase ion-molecule 
and neutral-neutr al chemistry {e.g. |Herbst & Leung |1989| 



Smith et al. 2004 1. Large abundances of carbon-chain-bearing 
molecules are often taken as a sign of chemical youth in 



molecular clouds (e.g. |Suzuki et al.|1992"{ Hirota et al.| 
Takakuwa et al. 201 ill. Indeed, the large abundances of these 
species in TMC-1 can be explained as being due to the 'early- 
time peak' that occurs when reactive carbon is still freely 
available in the gas (before it becomes almost completely 
locked up in CO). In environments denser than TMC-1, gas- 
phase species freeze out more rapidly onto dust due to in- 
creased collision rates, and this can become an important 
driver of large-molecule chemistry: Assuming equal sticking 
probabilities for all species involved in gas-grain collisions, 
lighter species (including atomic oxygen) travel at higher ve- 
locities and therefore freeze out more quickly than the heav- 
ier, slower-moving carbon chains. Because O is a primary 
destructive reactant for carbon-chain-bearing species, their 
abundances rise as oxy gen freezes out. This is kn own as the 
'freeze-out peak' {e.g. Brow n & Charnley]|1990| l, and may 
be a plausible explanation for the large carbon chain abun- 
dances observed in Cha-MMSl and other dense cores. Evi- 
dence in favour of a freeze-out-peak chemistry in Cha-MMS 1 
is given by the C O depletion and large NH3 abundance ( |Ten- 
nekes et al. 20 0~6]l. 



Alternatively, |Sakai et al.| (|2008b 2009a) hypothesised that 
elevated carbon chain abundances in L1527 arise as a result 
of heating by the newly-formed protostar IRAS 04368+2557. 
Once the outer envelope reaches a temperature > 30 K, it 
is theorised that grain-surface methane begins to sublimate, 
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Fig. 5. — Rotational diagrams for HC3N, HC5N, CCS and 0-0C3H2, including data from the present work and from Kontinen et al. ( 2000 1. Lines of best fit are 
shown (with corresponding LTE excitation temperatures T and column densities N). Integrated line intensities have been opacity-corrected (see Section[3.2.1(. 



which then reacts with gas-phase C + to form hydrocarbon 
ions. These subsequently engage in ion-molecule reactions 
and give rise to a so-called 'warm carbon-chain chemistry' 
(WCCC). Chemical models show that this mechanism is ca- 
pable of producing elevated abundances of unsaturated hydro- 
carbons including polyynes and cyanopolyy nes in the warm 
regions surrounding low-mass protostars ( |Hassel et al.| 2008 ). 
The detection of high-excitat ion-energy c arbon chain emis- 
sion lines in LI 527 by |Sakai et al.|p008b) is consistent with 
this theory. High-excitation lines have not yet been detected 
in Cha-MMS 1 , but we are presently undertaking new observa- 
tions to search for the m. The HC3N 7 = 10 - 9 map of Cha- 



MMS 1 published by Kontinen et al. (2000) shows a north- 
south elongation in the peak emission contour, which covers 
both the location of the protostar and the more southerly peak 
in our observed HC3N 7 = 4 — 3 map. The 7 = 4 — 3 emission 
originates from a rotational energy level (E u ) only 4.3 K above 
the ground state, whereas the 7 = 10-9 emission originates 
from E u — 24 K and therefore provides evidence for a signif- 
icant contribution to the HC3N emission from warmer, more 
excited gas heated by the protostar. It is possible that carbon 
chain abundances may b e enhanced as a re sult of WCCC in 

IpOOOl did not detect any 
32 



this region. However, |Kontinen et al. 
emission from the high-excitation {E u 



52 K) 7 = 15 - 14 



line of HC3N, which indicates that any warm and dense region 
must be small compared with the 37" SEST beam. 

According to the physical m odel of protostar evolution by 
Masunaga & Inutsuka (2000), during the FHSC stage soon 
before the protostar begins fusion, the temperature reaches 
above 30 K at radial distances from the core R30 < 20 AU. 
The protostar subsequently ignites and warms the surround- 
ing envelope over time. After 2 x 10 4 yr, is predicted to 
increase to ~ 1000 AU, and increases further during the mid- 
to-late stages of the Class phase, but by a relatively smaller 
amount. The Mopra beam HWHM of as 40" corresponds to 
6000 AU at the distance of Cha-MMS 1, so emission from the 
cool, extended outer envelope of the (collapsing) protostar 
dominates the observed spectra, and any emission from the 



compact methane-sublimation region (inside /?3o) would be 
severely diluted regardless of the evolutionary state of the pro- 
tostar, making the WCCC signature difficult to detect. In the 
case of Cha-MMS 1 (with its very low luminosity), the warm 
region close to the protostar is probably too small to make a 
significant contribution to our observed emission. Our non- 
detection of the more highly-excited 2 TLi/2 7 = 27/2 - 25/2 
emission lines of CeH (with E u = 31 K) is not surprising 
given the expected beam dilution. 

On the other hand, if Cha-MMS 1 contains a VeLLO un- 
dergoing episodic accretion, the size of the warm (methane- 
sublimation) region is predicted to vary in accordance with the 
accretion rate (see for example, Lee 2007 ). Therefore, if Cha- 
MMS 1 is presently at or near to a minimum in the accretion 
rate cycle, there may still be abundant carbon chains around 
the protostar, left over from a previous epoch of WCCC that 
occurred when the accretion rate and protostar luminosity 
were greater, and the size of R30 correspondingly larger. It 
must be emphasised, however, that the 'cold' carbon chain 
emission peaks (see Figures [2] to HI are offset from the proto- 
star centre, so WCCC is an unlikely explanation for the major- 
ity of the carbon-chain-rich material surrounding Cha-MMS 1 . 
Smaller telescope beam sizes combined with the observation 
of high excitation-energy lines to probe warmer gases will be 
required for further analysis of the carbon chemistry inside the 
putative methane sublimation zone and warm inner envelope 
of Cha-MMS 1. 

4.1.2. Anions 

The polyyne anions C 4 H~ and CeH were not detected. 
Three-sigma upper limits on their respective anion-to-neutral 
ratios are [C 4 H-]/[C 4 H] < 0.02% and [C 6 H-]/[C 6 H] < 10%, 
which include the errors in the observed abundances of the 
neutral species. The largest anion-to-neutral ratios observed 
in any astronomical source so far (for bo th species) have been 
in L1527 where [C 4 H']/[ C 4 H] = 0.01% jAgundez etal.|2008 1 
and [C 6 H-]/[C 6 H] = 9% ( |Sakai et al.|2008c) , which are con- 
sistent with our results. 
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In chemical models fo r dense clouds (e.g. Millar et al .|2007 i 
Harada & Herbst|2008 1, the observed C(,H~ anion-to-neutral 



ratios (on the order of a few percent) are reproduced with good 
accuracy. However, our non-detection of C4H provides a 
new example of the observed C4H anion-to-neutral ratio be- 
ing substantially less than predicted by theory (see Herbst & 
|Osamura|2008] ). 

4.1.3. Methanol 

Methanol is found in the gas phase in dense molecular 
clouds, prestellar cores and in the vicinit y of protostars, with 
column densities ~ 10 13 - 10 ls cm" 2 (e. g. Friberg et al. 1988; 
|Maret et al.1[2O05) |Buckle et aLp056) . However, most of 
the methanol in col d, den se protostellar envelopes is found 
in the form of ice ( |Boogert et al.||2008|l. Its formatio n via 
gas-phase chemistry is inefficient (Geppert et al. 2006}, and 
methanol is widely believed to be formed by successive hy- 
drogenation of CO on cold dust grain surfaces. Its appearance 
in the gas phase in cold regions has been hypothesised to be 
due to cosmic -ray-induced desorption, reactive desorption, or 
heating in grain-grain collisions, whereas inside protostellar 
envelopes, heating of the material to temperatures > 100 K 



results in complete evaporation of the ices (e.g. Rodgers & 
Charnley 2003 ). 

The CH3OH map in Figure |3] does not rule out the exis- 
tence of methanol in warm gas in the vicinity of the protostar, 
but similar to the carbon-chain-bearing species, the fact that 
the emission peak is offset by ~ 10 4 AU shows that protostel- 
lar heating probably does not explain the majority of the ob- 
served methanol. The lack of spatial correlation between the 
CH3OH and carbon chain distributions matches the observa- 
tions made by Buckle et al. (2006 ), who used chemical models 
to infer the physical and chemical histories of HC^N-rich and 
CH30H-rich gas clumps. The observed methanol peak off- 
set in Cha-MMSl is co nsis tent with the theory presented by 
IMarkwick et al.| ( |2000| > and |Buckle et al.| ( |2006) , that Alfven 
waves emitted by nearby collapsing protostars cause transient 
removal of ice mantles from the surrounding dust, resulting in 
the liberation of methanol into the gas phase at a distance from 
the protostellar core on the order of th e MHD-wave damping 
length (~ 10 4 AU). |Markwick et al.| ( |2000| used this theory 
to successfully explain the origin of the abundance gradients 
of various species observed along the ridge of TMC-1. The 
young age of Cha-MMSl implies that protostellar collapse 
is currently in progress, resulting in the production of Alfven 
waves that could be responsible for desorption of the observed 
methanol. 

From the ratios of our observed CH3OH lines with those of 



Kontinen et al. (2000 ), differing temperatures were calculated 
for the A and E methanol nuclear spin species of 7.1 ± 0.9 K 
and 4.7 + 0.6 K, respectively. The column densities of A and 
E methanol (given in Table [2} also differ, but within the error 
bars. The ratio of the sums over the statistical weights of the 
quantum states of A and E methanol approaches unity, so in 
the limit of high temperature (in thermodynamic and chemical 
equilibrium), the abundances of A and E methanol should be 
approximately equal. Convergence of the A and E abundances 
is calculated to occur for t emperatures above about 30 K (see 
|Wirstr6m||2009| |Wirstrom et al.||201l| . At lower tempera- 
tures, A-methanol is the more abundant species because its 
ground state is 7.9 K lower in energy than that of E-methanol. 
Given that both spin species are likely to have the same spa- 
tial distribution in the ISM, their relative column densities in 
Cha-MMS 1 may be more accurately compared by assuming a 



common excitation temperature. In LTE at 6.1 K, the column 
densities of A and E methanol calculated from our observa- 
tions are (7.2 + 0.7) x 10 12 cirT 2 and (6.6 ± 0.7) x 1 12 cm" 2 
Using only t he higher- frequency lines observed by |Kontinen| 
etaLl (2000), we calculate the respective column densities 



to~be (73± 0.7) x 10 12 cm" 2 and (7.9 ± 0.8) x 10 12 cm 
These results are consistent with equal abundances of A and 
E methanol, implying formation at close to statistical equilib- 
rium. 



4. 1 .4. Deuterium fractionation 

In cold interstellar clouds, molecular ions become enriched 
in deuterium as a result o f exotherm ic gas-phase fractionation 
re actions involvin g HD (Milla ret al.|1989| ). 

|Belloche et al.| ( |2006| ) found a deuterated N 2 H + fraction 
in Cha-MMSl of [N 2 D + ]/[N 2 H + ] = 11%, which is typi- 
cal for prestellar cores and young protostars observed in 
the nearby Galaxy. The cyanoacetylene deuteration fraction 
[DC 3 N]/[HC 3 N] = 3.6 ± 1.5% is within the range of values 
previously observed in TMC-1 (1 - 6%; |Turner 2001] [SaTto| 



et al. 



2002 1, and is similar to the carbon-chain-rich prestellar 



core L1544 (6%; Howe e t alJ1994) , and the protostar L1527 
(3%; [Sakai et al. 2009b). These results are con sistent with 
chemical models for cold interstellar clouds (e.g. |Roberts & 
Millar|2000||Turner|2001HPagani et al.|2009] ). 



4.1.5. Cyclopropenylidene 

A large 0C3HD column density of 4.5^p 3 x 10 12 cm" 2 is 
found in Cha-MMSl, which is about 5 times greater than 
in LI 527 and an order of magnitude greater than in TMC-1 
( Sakai et al.|2009b Turner|2001[ ). Accounting for the rather 
large uncertainty in the C-C3H2 column density (see Section 
|3.2.3| ), we calculate that the cyclopropenylidene deuteration 
fraction [c-C 3 HD]/[c-C3H 2 ] is between 5% and 34%. This is 
not significantly different from the range of values (5 - 15%) 



reporte d in a sample of twelve dark clouds by |Bell et al. 
( 1988 ), and the value of 7% reported in L1527 by Sakai et al. 



(2009b), but is consistent with previous results which show 
the c-C3H 2 deuteration fraction is larger than that of HC3N 
(e.g. ITurner 2001 ; Sakai et al. 2009b). Chemical models pre- 
dict comparable levels of deuteration for 0C3HD and HC3N; 
deuterium enrichment in both these species is consi dered to 
originate predominantly from the CH 2 D + ion (e.g. Roberts 
|& Mil lar 2000; Turner 2 001| ). Deuteration may also occur 
by deuteron transfer as a result of collisions with H 2 D + , fol- 
lowed by dehydrogenation upon recombination. Discrepancy 
between the deuteration fractions of 0C3HD and HC3N is 
indicative of an incomplete understanding of the deuterium 
chemistry of these two species. The fact that 0C3H2 can be- 
come doubly deuterated may partly explain why 0C3HD is 
so abundant: if sufficient c-C3D 2 forms, it would act as a 
buffer for 0C3HD in subsequent proton-deuteron exchange 
processes. It may thus be of interest to search for c-C3D 2 in 
Cha-MMSl and similar environments. However, we are not 
aware of any published laboratory studies on the rotational 
spectrum of c-C3D 2 ; such a study would be warranted. 

The C-C3H2 molecule exists in two nuclear spin configu- 
rations: ortho and para, and it has been suggested that the 
ortho-to-para ratio (OPR) may be us ed as a probe of interstel- 
lar cloud ages (Morisawa et al.|2006[ ). However, our present 
data are insufficient to usefully constrain the OPR. Additional 
observations of c-C3H 2 and 0C3HD lines from a range of en- 
ergy energy levels and with a range of optical depths will be 
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required in order to better constrain the column densities of 
these species and derive more reliable values for the 0C3H2 
OPR and deuteration fraction. 

4.1.6. Sulphuretted species 

We observed large abundances of the sulphur-bearing car- 
bon chains C„S (« = 1-3), on the same order of magnitude as 
found in various parts of the Taurus molecular cloud complex 
(e.g. [Hirota et al.|2004| >. Our [CCS]/[C 3 S] column den sity ra- 
tio (5 .7 + 2.1) is greater than the values observed by |Hirota| 
[etaLl of 2.1, 2.3, 2.0 in the dense cores L1495B, L1521B and 
L1521E, respectively, and is more similar to the value of 5.1 
found in TMC-1. The relatively large value in Cha-MMSl is 
primarily attributable to a relatively low C3S column density, 
which is more consistent with the lower C3S column densi- 
ties observed in the carbon-chain-rich cloud cores L1512 and 



LI 251 A by Cordine r et al | ( pOTT i. 

Sulphur is believed to be incorporated into gas-phase hy- 
drocarbon molecul es in the ISM mainly b y ion-neutral chem- 
istry involving S + ( |Millar & Herbst|1990| l. Recombination of 
HC„S + species then produces sulphuretted carbon chains. The 
main source of HC3S + is from the reaction S + + 0C3H2 — > 
HC3S + + H, followed by recombination of HC3S + to produce 
C 3 S. Thus, the relatively low C 3 S abundance in Cha-MMSl 
could be explained by a relative underabundance of 0C3H2 
compared with L1495B, L1521B and L1521E. An alternative 
method for the synthesis of C3S is by the neutr al-neutral reac- 
tion C 2 H + CS — > C3S + H ( |Smithetal.|2004l l. However, this 
reaction is relatively slow (two orders of magnitude slower 
than the ion-neutral reaction), and is therefore expected to be 
less important in the production of C3S. 

The [N2H + ]/[CCS] ratio may be used as a measure of chem- 
ical age in interstellar clouds on the basis that carbon chains 
are most abundant at early times, whereas N2HL becomes 
more abun dant later on as CO freezes out ( e.g. |Tak akuwa 
et al.|201 7) . Using the N 2 H + abundance from Belloche et al. 
(2006), we derive [N 2 H + ]/[CCS] = 0.5 ± 0.3 in Cha-MMSl, 
which is similar to the value of 0.52 ±0. 07 found in the y oung, 
carbon-chain-rich VeLLO L1521F (Takakuwa et al. 201 1 ), in- 
dicating a comparable chemical age, but less than the value 
1.18 + 0.21 in the more evolved VeLLO IRAM 04191. 

Only a lower limit could be placed on the CS column den- 
sity because the optical depth of the 7=1-0 lines approaches 
infinity for the adopted temperature range of 6.1 ± 1.4 K. 

4.1.7. Other molecules 

Acetaldehyde (CH3CHO) was detected in Cha-MMSl 
with a column density of 3.6 x 10 12 cm 2 . It is an im- 
portant molecule in terrestrial organic chemistry, and has 
been detected in star-forming regions, protostars and translu- 
cent/dense molecular clouds. A range of gas and ice-phase re- 
action pathways have been proposed to explain its abundance 
in space (see [Charnley 2004 1 Benn ett et al.||2005] and refer- 
ences therein). However, the size of the Mopira beam does not 
place any useful constraint on the origin of the emission rel- 
ative to the protostar, and therefore, apart from highlighting 
the diversity of the organic chemistry of this region, the ob- 
served CH3CHO abundance of ss 10~ ll) does not significantly 
constrain the physics or chemistry of the molecule or the pro- 
tostar. Our observations are, however, consistent with equal 
abundances of gas-phase A- and E- CH3CHO, indicating sta- 
tistical equilibrium. 

Ethanol (CH3CH2OH) has previously been observed with 
relatively large abundances of ~ 10~ 8 in 'hot core' regions 



surrounding high-mass protostars (e.g. [Millar et al. 1995), 
but is not present in quiescent (non-star-forming) molecular 
clouds. There have been two reported detections to-date in 
low-mass protostars: an abundance of ~ 10 ~ 8 was found in 
the hot corino IRAS 16293-242 2A bv|Huang etlI1 ( |2005l l and 
~ 10" 7 in the LI 157 outflow by |Arce etaLl (2008) ' Wederive 
a (beam-averaged) upper limit on the ethanol abundance of 
5xl0- I0 in Cha-MMSl. 

We also searched for the protonated carbon dioxide ion 
HCOj and derive an upper abundance limit of 4.7 x 10" 11 . 
This molecule was detected by |Sakai et al.| ((2008a) in the 
vicinity of the low-mass protostar L1527, with a beam- 
averaged abundance of 1.3 x 10 -12 . HCOJ is believed to be 
formed predominantly by protonation of CO2, and the gas- 
phase CO2 abundan ce can be estimated using Equation 5 of 
Sakai et al.||2008a| l. Using a CO abundance of 5 x 10~ 5 (cal- 
culated using t he C 18 Q column density of 2.5 X 10 15 cm' 2 
reported by |Kontinen et al.|2000| and an intrinsic 16 0/ 18 ra- 
tio of 500), we derive an upper limit of n(C02) < 4 x 1 0~ 6 «h,, 
averaged o ver the Mopra beam. Spitzer observations (Furlan 
|et al.|20 08) show a solid CO2 abundance of 5 x 10~ 6 towards 
L1527 and a similar ice abundance may be expected in the 
cool, outer envelope of Cha-MMS 1 . The relatively low gas- 
phase CO2 in Cha-MMS 1 implied by our observations, com- 
bined with the upper limit on the ethanol abundance, shows 
that the energetic processing of matter that would cause the 
sublimation of dust grain mantles in the protostellar envelope, 
due to heating by the protostar or an outflow shock, for exam- 
ple, must be lacking in Cha-MMSl, or at least confined to a 
relatively small area compared with the « 90" Mopra beam. 

4.2. Comparison with other sources 

The molecular abundances in Cha-MMSl are indicative 
of an active (non-equilibrium) carbon chemistry, resulting in 
the synthesis of large quantities of unsaturated carbon-chain- 
bearing species such as the polyynes and cyanopolyynes. The 
carbon-chain abundances are similar to those in young, qui- 
escent cores such as T MC-1 and the 'carbon-chai n-producing 
regions' obse rved by Suzuki et al.| (|1992|l and |Hirota et al. 



(2004, 2009). However, as pointed out by Tennekes et al. 
(2006), the CO depletion and large NH3 abundance in Cha- 
MMS1 are hallmarks of a later stage of chemical evolution, 
associated with the freezing-out of gas-phase species onto the 
dust, as has been observed in prestellar cores and the outer en- 
velopes of other young, low-to-inte rmediate-mass protostars 
flCaselli et al.|1999||j 0rgensen et al |2004| |Crapsi et al.|2005| 
Alonso-Albi et al.|201o| ~ 

Only a few detections of long carbon-chain-bearing species 
in low-mass protostellar envelopes have been reported to-date. 
The smaller species HC 3 N, C 3 H 2 and C 4 H in Cha-MMSl are 
a factor of a few to an order of magnitude more abundant than 
those surro unding the VeLLO LI 52 IF and the prestellar core 
L1544 (seejSuzuki et al |[T992j |Gupta et al.||2009l |Takakuwa| 



et al. 2011 1, both of which are considered to be relatively 
carbon-chain-rich sources. Due to their evolved state, the ap- 
pearance of carbon-chain-bearing molecules in these sources 
is likel y to be the result of 'freeze-out-peak' chemistry (see 
Section |4TTTj ). 

Carbon chains have also been hypothesised to arise in large 
abundances in protostellar envelo pes as a result o f warm 
carbon-chain chemistry (WCCC) ( |Sakai et al.|2009a) . Of aU 
the dense molecular clouds in the literature for which multi- 
ple complex molecules have been observed, the abundances in 
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Cha-MMS 1 match most closely with those in the well-studied 
Class 0/1 protostar L1527 in the Taurus molecular cloud com- 
plex ( [Sakai et al.||2008b[ |2009a|b| l, where WCCC has been 
theorised to be occurring. Indeed, the abundances of all ob- 
served species in Cha-MMS 1 are within an order of magni- 
tude of those observed in L1527. 

The gas densities in the envelopes of Cha-MMS 1 and 
L1527 are quite similar (~ 10 6 cm -3 ), so the similarities in 
their chemistries may be related to similarities in the physical 
conditions. It is tempting, therefore, to hypothesise that the 
large abundances of organic molecules around Cha-MMS 1 
also arise as a result of WCCC. However, Cha-MMS 1 is likely 
to be at a significantly earlier stage of evolution than L1527, 
and has an estimated luminosity ~ 1/100 as large (J0rgensen 
|et al.|2002||Belloche et al.|201 1|>. The low IR flux and lack of 
outflow fro m the protostar ( Belloc he et al.|[2006{ |Hira matsu 
|et al.|2007| indicate that Cha-MMS 1 is very close to the be- 
ginning of the Class phase, whereas L1527 is more evolved, 
with bright scattered IR emission and clear outflow cavities, 
consistent with an edge-on C lass I protostar obscured by its 
own disk (|Tobin et al.|20Q8] >. Assuming a uniform accretion 
rate, |Ohashi et al.| ( |199 7[ l calculated the age of the L1527 pro- 
tostar to be ~ 10 5 yr, which would place it at the end of the 



Class ph ase (based the Class lifetime derived by Enoch 
|et al.|2009] >. Chemical timescales in dense interstellar clouds 
(in particul ar, molecu lar freeze-out timescales; see Equation 4 
of |Buckle et al.|2006| are significantly less than the estimated 
-10^ yr age difference of these two protostars. Therefore, the 
chemical similarity between Cha-MMS 1 and LI 527 may be 
coincidental. Notwithstanding the possibility of episodic out- 
bursts, the low luminosity of Cha-MMS 1 implies a very small 
hot-core region (probably < 10 AU in size), that would not 
be detectable with single-dish microwave telescopes. Instead, 
the large carbon chain abundances observed in the envelope 
of Cha-MMS 1 may be present because this young protostar is 
formi ng from a very chemically-rich parent molecular cloud. 

Crapsi et al. p005) derived CO depletion factors Sco ~ 10 
in L15 21F and L1544. Using the CO abundance from Section 
|4.1.7| we calculate <5co ~ 2 in Cha-MMS l (assuming an un - 



Crapsi et al.|2005| >. 



depleted CO abundance of 9.5 x 10~ 5 after|( 

The size of the warm region around the Cha-MMSl proto- 
star within which CO ice is desorbed is likely to be negligi- 
ble compared with the telescope beam, so low CO depletion 
may be a sign that the freezing-out of gas-phase species onto 
the dust is less far advanced. This would be consistent with 
a relatively short prestellar phase for Cha-MMSl. Alterna- 
tively, the low CO depletion could be a consequence of grain 
mantle destruction; significant turbulent motion of the gas 
around Cha-MMS 1 may be inferred from the relatively broad 
molecular linewidths (« 0.78 kms~'), which are more than 
an order of magnitude greater than the thermal limit in most 
cases, and are s ignificantly broade r than observed in L1521F 



(< 0.3 kms" 1 ; Gupta et al. 



2009 1. This may be a signature 



of energetic cloud motions that cause collisions between dust 
grains, resulting in the liberation of mantle species into the gas 
phase. Possible physical processes responsible for mantle de- 
struction include inter-clump collisions or Al fven w aves emit- 
ted by the protostar. As discussed in Section [4.1.3| the libera- 
tion of mantle species ma y explain the origin of the observed 
methanol. According to Markwick et al. (2000), grain man- 
tle desorption also gives rise to large gas-phase abundances of 
polyynes and cyanopolyynes. 
To help ascertain the origin of complex molecules in the 



vicinity of Cha-MMSl, direct observations of the core region 
at higher spatial resolution are required, which will also help 
to establish the precise physical nature of the source and iden- 
tify molecular signatures of the earliest stages of star forma- 
tion. For example, observations of the gas-phase products of 
grain-surface chemistry (in particular, alcohols, saturated or- 
ganic molecules and deuterated sp ecies) will probe th e pres- 
ence of a possible hot corino (e.g. Cazaux et al.|2003 1. High- 
resolution maps will permit the search for molecular abun- 
dance gradie nts due to Alfven waves predicted by |Markwick 
|et al.| l |2000) . The identification of concentrated abundances 
of carbon chains around the protost ar would provide ev id ence 
for wa rm carbon chain chemistry (Sakai et al. 2010). Lee 
(2007) theorised that by mapping the JSl2H + and CO abun- 
dances at high angular resolution, and using infrared CO2 ice 
observations in the protostellar envelope, it should be possible 
to observe the chemical effects of the episodic ac cretion pro- 
cess hypothesised to be present in VeLLOs. |Omukai| (2007) 
calculated that high-excitation submillimeter H2O emission 
lines may be used as a tracer for the presence of the infall 
shock-front that is predicted to be a characteristic property 
of first hydrostatic cores. A compact, low-velocity outflow 
is predicted to be another cha racteristic of the FHSC phase 
( [Machida et al.||2008| |Tomida et ail2010| >. It may thus be 
possible to determine whether Cha-MMSl is a VeLLO or an 
FHSC using high angular-resolution (~ 0.1 - 1") molecular 
line observations. 

5. CONCLUSION 

The envelope of the low-luminosity protostar Cha-MMSl 
exhibits an unusually rich chemical diversity. It contains very 
large abundances of organic molecules, indicative of active 
gas-phase and grain-surface chemical processes near to the 
protostar. The abundances of long carbon chains (including 
CgH and HC7N) are particularly large, and similar to those 
found in the most carbon-chain-rich regions of the Galaxy. 
These observations show that complex organic molecules and 
carbon-chain-bearing species can be abundant in star-forming 
gas, and may therefore be of relevance to studies of the chem- 
ical evolution of protoplanetary disks. The molecular abun- 
dances in Cha-MMSl are similar to those in the envelope 
of the older, more luminous Class 0/1 protostar L1527, and 
are greater than have been reported so far for any other low- 
luminosity, low-mass Class protostellar envelope (includ- 
ing the VeLLOs L1521F, L1014 and IRAM 04191). This 
may be attributable to a combination of factors including (1) 
a young chemical age for the Cha-MMSl natal interstellar 
cloud, (2) freezing out of atomic oxygen onto the dust, and (3) 
the liberation of grain mantle species into the gas by various 
mechanisms such as clump-clump collisions, Alfven waves or 
warming by the protostar. 

The observations presented here are of lines from relatively 
low excitation levels which trace predominantly cool mate- 
rial; the chemistry of the warmer regions closer to the proto- 
star may be examined by the observation of higher-excitation 
lines at higher spatial resolution. This will permit further anal- 
ysis of the relationship between protostar and envelope chem- 
istry. Establishing the molecular inventory of the envelope 
of Cha-MMSl and other young protostars will be of funda- 
mental importance towards the goal of understanding the ini- 
tial chemical reagents available for the formation of planetary 
systems. 
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